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ABSTRACT 

We have undertaken an extensive study of X-ray data from the accreting millisecond pulsar 
XTE J175 1-305 observed by RXTE and XMM-Newton during its 2002 outburst. In all aspects 
this source is similar to a prototypical millisecond pulsar SAX J1808. 4-3658, except for the 
higher peak luminosity of 13 per cent of Eddington, and the optical depth of the hard X-ray 
source larger by factor ^2. Its broad-band X-ray spectrum can be modelled by three compo- 
nents. We interpret the two soft components as thermal emission from a colder (kT ~ 0.6 
keV) accretion disc and a hotter (~ 1 keV) spot on the neutron star surface. We interpret the 
hard component as thermal Comptonization in plasma of temperature ~40 keV and optical 
depth of ^1.5 in a slab geometry. The plasma is heated by the accretion shock as the material 
collimated by the magnetic field impacts on to the surface. The seed photons for Comptoniza- 
tion are provided by the hotspot, not by the disc. The Compton reflection is weak and the disc 
is probably truncated into an optically thin flow above the magnetospheric radius. Rotation of 
the emission region with the star creates an almost sinusoidal pulse profile with rms ampli- 
tude of 3.3 per cent. The energy-dependent soft phase lags can be modelled by two pulsating 
components shifted in phase, which is naturally explained by a different character of emission 
of the optically thick spot and optically thin shock combined with the action of the Doppler 
boosting. The observed variability amplitude constrains the hotspot to lie within 3^1° of the 
rotational pole. We estimate the inner radius of the optically thick accreting disc of about 40 
km. In that case, the absence of the emission from the antipodal spot, which can be blocked 
by the accretion disc, gives the inclination of the system to be > 70°. 
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1 INTRODUCTION 

Millisecond radio pulsars are believed to be spun-up to their fast ro- 
tational speeds during the p hase of accretion fr om their companion 
stars (for a review see e.g. Bhattachar vall 1 995ft . Though fast radio 
pulsars have been known for over 20 years, accretion-powered mil- 
lisecond pulsars had proven to be an elusive missing link for a long 
time. SAX J1808. 4-3658, the first millisecond pulsar (spin period 
P s = 2.5 ms) in a low-mass X-ray binary (LXMB) was discovered 
in 1998 (Wiinands & van der Klis 1998a) by Rossi X-ray Timing 
Exp lorer (RXTE). Recent di scoveries of XTE J1751-305 (P a = 2.3 
ms: IMarkwardt et alJl2002l here after M02), XTE J0902-314 (P s 
= 5.4 ms:lGallowavetalJl2002h. XTE J1807-294 (P s = 5.3 ms; 
Markwardt. Smith & Swankl2003l) . XTE J1814-338 (P s = 3.2 ms; 
Mark wardt & Swankl2003l) . and IGR J00291+5934 (P a = 1.67 ms; 
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lEckert et al1l2004t IMarkwardt et alJl2004l) have brought the total 
number of currently known accreting millisecond pulsars to six. 

iGierlinski. Done & Barrel feOOl henceforth GDB02) anal- 
ysed energy spectra of SAX J1808.4-3658 from the 1998 out- 
burst. The X-ray spectral shape remained almost constant through- 
out the outburst and could be fitted by a blackbody and thermal 
Comptonization. Weak Compton reflection (Q/2n ~ 0.1) was also 
detected. GDB02 interpreted the blackbody as emission from the 
heated spot on the neutron star, and the Comptonization taking 
place in an accretion column, as the material collimated by the 
magnetic field impacts on to the neutron star surface. The phase- 
resolved X-ray spectra showed that pulsatio n of these t wo compo- 
nents is shifted in phase. Poutan en & Gierlihskil 12003. henceforth 
PG03) applied the full relativistic treatment to the energy-resolved 
pulse profiles of SAX J1808.4-3658 and studied the effects of dif- 
ferent angular distribution of the emission pattern from the black- 
body and Comptonization components. They estimated the neutron 
star radius to be ~ 6.5 and ~ 1 1 km for a 1.2 and 1.6 M© neutron 
star, respectively. 

XTE J1751-305 is in a very tight binary, with the orbital 
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period of only 42 minutes and the mass function of the pulsar 
1.3 x 10 -6 M©, which gives a minimum mass for the compan- 
ion of about 0.014 Mq (M02). The distance, D, to the source is not 
known, though as it might be related to the Galactic centre, there- 
fore we assume D ~ 8.5 kpc. M02 estimated the distance to be 
greater than 7 kpc, using indirect arguments. 

In this paper we examine the X-ray spectra of XTE J175 1-305 
from its April 2002 outburst, observed by XMM-Newton and RXTE. 
In Section[4]we study the (phase-averaged) energy spectra, in par- 
ticular the broad-band 0.7-200 keV simultaneous spectrum from 
both instruments. We model them by the blackbody and Comp- 
tonization from hotspot and accretion shock, respectively. We also 
show that another soft component, from the accretion disc, is re- 
quired. In Sections|5|and|6|we study phase-resolved spectra and in- 
terpret them in terms of two independently pulsating spectral com- 
ponents. In Section|8|we discuss our results and derive constraints 
on the accretion flow geometry which occur to be very similar to 
those found for SAX J1808.4-3658 (GDB02). 



2 OBSERVATIONS 
2.1 RXTE 

RXTE observed XTE J175 1-305 from 3 rd to 30 th of April 2002. 
We reduced these data using FTOOLS version 5.3. Proportional 
Counter Array (PCA) units 1 and 4 were switched off for most 
of this period. As unit lost its propane layer in May 2000 and 
its calibration is uncertain, we decided to use the data from units 2 
and 3 only. For each of the observations with the unique identifier 
(obsid) we extracted energy spectra in 3-20 keV band. Analogous 
Crab spectra from the same period show features at a level of 1 per 
cent, most likely of the instrumental origin. Therefore, we applied 
1 per cent systematic errors in each of the PCA energy channels in 
addition to statistical errors. 

In addition to the PCA data we also extracted 20-200 keV 
High-Energy X-ray Timing Experiment (HEXTE) spectra from 
both detector clusters. During observation of April 3 HEXTE did 
not perform 'rocking' and did not record background, so we could 
not use these data. Table \\\ contains log of the RXTE observations 
analysed in this paper. 

XTE J175 1-305 is only 2° away from the Galactic centre, 
so its X-ray spectrum is contaminated by the diffuse Galactic ridge 
emission. Fortunately, RXTE performed several observations in the 
direction of XTE J1751-305 after the outburst (see fig. 1 in M02). 
Between 20 th and 25 th of April PCA count rate remained constant 
and minimal (there was a small outburst later on), indicating its 
origin from the background. We accumulated average PCA spec- 
trum from this period (with average count rate of 12.01 ± 0.04 
s _1 in detectors 2+3, all layers) and used it as an additional back- 
ground for outburst observations. TableQcontains PCA count rates 
corrected for this effect. HEXTE spectra were not affected by the 
diffuse emission. 

We extracted phase-resolved energy spectra from the PCA 
Event mode files (configuration E_125us_64M_0_ls) with timing 
resolution of 122 fis (except for the April 3 observation where we 
used GoodXenon mode files, with resolution of 1 ps). We gener- 
ated folded light-curves in 16 phase bins, for each PCA channel, 
using events from all layers of units 2 and 3. We have chosen be- 
ginning of the phase (<j> = 0) at the bin with lowest 3-20 keV count 
rate. Photon arrival times were corrected for orbital movements of 
the pulsar and the spacecraft, using ephemeris from M02. Back- 
ground files were created from standard models, with addition of 
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Figure 1. 3-20 keV unabsorbed flux of observations 1-21 (TableQJ. The 
solid lines correspond to F oc e~ */ 7 ' 2 and F oc e~ '/0-63 . 

the April 20-25 spectrum to represent the diffused emission. We 
also extracted power density spectra (PDS) from the same event 
data files, but using units 0, 2 and 3. We have created power density 
spectra in the 1 /i28-512 Hz frequency range from averaging fast 
Fourier transforms over 128-s data intervals. 

2.2 XMM-Newton 

XMM-Newton observed XTE J1751-305 on 2002 April 7, with 
the total exposure of 33 ks. We reduced the data from Euro- 
pean Photo Imaging Camera PN (EPIC-pn) detector using SAS 
version 5.4.1. and following guidelines in the MPE cookbook 
(http://wave.xray.mpe.mpg.de/xmm/cookbook I. The detector oper- 
ated in timing mode where spatial information was compressed into 
one dimension. We extracted EPIC-pn spectrum from a stripe in 
raw CCD coordinates RAWX (36,56), collecting single and double 
events only. The background was extracted from the adjacent stripe 
in RAWX (16,36). 

All spectral analysis (both phase-resolved and phase- 
averag ed) wa s done using the XSPEC 11.2 spectral package 
I Arnaud 1996). The error of each model parameter is given for a 90 
per cent confidence interval, except for the pulse profiles in Sec. [6] 
where we used 1<t errors. The relative normalization of the PCA 
and HEXTE instruments is uncertain, so we allowed this to be an 
addition free parameter in all spectral fits. 



3 OUTBURST 

Fig.Qshows the evolution of the 3-20 keV unabsorbed flux during 
the outburst (see also M02). The peak flux during observation 2 
was 1.5 x 10 -9 erg s _1 cm~ 2 . Using a Comptonization model 
(see Section[4j, we estimated the peak bolometric X-ray/7-ray flux 
of 3.1 x 10~ 9 erg s _1 cm~ 2 . For a distance of 8.5 kpc (assuming 
the source to be close to the Galactic centre) this corresponds to 
the bolometric luminosity of 2.7 x 10 37 erg s _1 , or 13 per cent 
of Eddington luminosity LEdd for a 1.4 Mq neutron star. This is 
much brighter than SAX J1808. 4-3658 with the estimated peak 
bolometric luminosity of 2.2 per cent of LEdd (GDB02). 

The light curve of XTE J1751-305 shows a striking similarity 
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Table 1. Log of RXTE observations analysed in this paper. Start and end times are in the UT days of April 2002. Exposures are in seconds and (background 
subtracted) count rates are in counts per second from detectors 2 and 3 of PCA and separately from cluster and 1 of HEXTE. During the first observation 
HEXTE did not observe background. 



PCA (2+3) HEXTE HEXTE 1 

No. Obsid Start End Exposure Count rate Exposure Count rate Exposure Count rate 



1 


70134-03-01-00 


3.653 


3.720 


3168 


250.5±0.6 


N/A 


N/A 


2 


70131-01-01-00 


4.644 


4.912 


10960 


237.3±0.5 


3520 


18.7±0.2 


3461 


15.2±0.2 


3 


70134-03-02-00 


5.162 


5.179 


720 


226.5±0.8 


257 


19.8±0.9 


263 


14.4±0.8 


4 


70131-01-02-00 


5.535 


5.773 


6384 


214.8±0.5 


2039 


17.3±0.3 


2012 


14.0±0.2 


g 


701 31 -01 -03-01 


6 525 


6 542 


1456 


i on f)+0 6 


425 


16.4±0.6 


431 


12.3±0.5 


6 


70131-01-03-00 


6.591 


6.763 


9200 


185.9±0.4 


2859 


15.3±0.2 


2791 


12.0±0.2 


7 


70131-01-04-00 


7.515 


7.753 


12416 


163.3±0.4 


3839 


13.7±0.2 


3809 


10.7±0.2 


8 


70131-01-05-03 


8.061 


8.083 


1168 


162.4±0.6 


432 


14.3±0.7 


410 


8.9±0.6 


9 


70131-01-05-02 


8.131 


8.149 


864 


154.8±0.7 


313 


14.4±0.9 


302 


10.3±0.8 


10 


70131-01-05-01 


8.199 


8.215 


576 


155.4±0.8 


204 


12.6±1.1 


198 


9.3±1.0 


11 


70131-01-05-04 


8.267 


8.467 


2368 


147.7±0.5 


754 


13.8±0.5 


749 


9.4±0.4 


12 


70131-01-05-000 


8.505 


8.743 


12544 


144.9±0.4 


3897 


13.4±0.2 


3799 


9.6±0.2 


13 


70131-01-05-00 


8.772 


9.007 


4752 


143.1±0.4 


1645 


12.9±0.3 


1586 


8.9±0.2 


14 


70131-01-06-00 


9.435 


9.601 


8800 


129.0±0.3 


2831 


12.6±0.2 


2731 


8.6±0.2 


15 


70131-01-06-01 


9.627 


9.799 


8464 


126.1±0.3 


2714 


12.1±0.2 


2660 


8.0±0.2 


16 


70131-01-07-00 


10.418 


10.525 


6592 


112.1±Q.3 


2093 


10.4±0.3 


2092 


7.1±0.2 


17 


70131-01-07-01 


10.551 


10.789 


11552 


108.6±0.3 


3619 


9.9±0.2 


3573 


6.8±0.2 


18 


70131-01-08-000 


11.409 


11.713 


15552 


94.3±0.3 


4840 


8.9±0.2 


4794 


5.9±0.1 


19 


70131-01-08-00 


11.742 


11.912 


6704 


91.2±0.3 


2376 


8.7±0.2 


2334 


5.5±0.2 


20 


70131-01-09-000 


12.333 


12.638 


9632 


56.7±0.2 


2934 


5.9±0.2 


2870 


3.3±0.2 


21 


70131-01-10-00G 


13.322 


13.496 


5120 


13.7±0.2 


1683 


1.3±0.3 


1686 


0.6±0.2 



to that of SAX J 1808. 4-365 8 (see fig. 1 in lGilfanov et alj 1998k see 
also M02). After the peak the flux declined exponentially, with e- 
folding factor of 7.2 days (10 days in SAX J1808.4-3658), until it 
reached a break after which the flux dropped suddenly with the e- 
folding factor of ~0.6 days (1.3 days in SAX J1808.4-3658). The 
break in the light curve is most likely associated with the onset of 
the cooling wave in the accretion disc (GDB02). 



4 PHASE-AVERAGED ENERGY SPECTRA 
4.1 Physical picture and spectral models 

If the neutron star magnetic field is strong enough, the accreting 
material will follow magnetic field lines and form a shock close 
to the neutron star surface (for the ac cretion geometry see fig. 12 
of GDB02 and Roma nova et all20 04). The shock is pinned down 
to the stellar surface at luminosity of a few per cen t of ^Edd 
iBasko & Sunvaevl Il976t iLvubarskii & Sunvaevl 1 19821) . Gravita- 
tional energy is dissipated in the shock and is transferred from 
the protons to the electrons by Coulomb collisions. The main cool- 
ing mechanism is Comptonization of soft photons provided by the 
neutron star. The stellar surface is heated under the shock, while 
the hard X-rays, produced in the shock, can also irradiate the sur- 
rounding surface, so that the black body emission region can cover 
a somewhat larger area. 

Such a physical picture is in agreement with observa- 
tions of SAX J1808.4-3658 (GDB02, PG03). We presume that 
XTE J175 1-305 is a similar source to SAX J1808.4-3658. Thus 
we expect emission from the hotspot on the neutron star surface, 
which we model as a blackbody. We also expect emission from the 
accretion shock, which we model by thermal Comptonization. We 
consider two possible geometries of the shock and hotspot, shown 



in Fig. [2] For the phase-averaged fits we assume that the angle be- 
tween the rotation axis and the magnetic pole is small, so that the 
slab model of the shocked region viewed at a fixed angle is justified. 
We will show later in Sec. 18. 3.41 that this assumption is feasible. In 
addition to the above components we might expect thermal emis- 
sion from the accretion disc and Compton reflection from the disc 
and/or neutron star surface. 

For the spectral description we use the models summarized 
in Table|2| Their spectral components are as follows. BBODYRAD 
is a single-temperature blackbody, normalized to its apparent area 
(at a given distance), which represents the emission from a hotspot 
at a neutron star surface. DISKBB is a multicolour (accretion) disc 
model iMits uda etalll984l) . normalized to the inner disc radius (at 
a given distance). We correct the inn er disc radius for the spec- 
tral hardening with factor / co i = 1.8 1 Shi mura & Takaharall 19951) 
and for the torque-free inner boundary condition, with factor £ = 
0.37 jGierlinski et ali n^^^: R ln = / c 2 iC-Rin, apparent- The exact 
value of the hardening factor is not well known, but recent nu- 
merical accretion disc simulations place it between 1.8 and 2.0, 
increasing with luminosity (Davis et al. 2004). Additional uncer- 
tain ty comes from general a nd special relativistic corrections (see 
e.g. ICunningham1ll975t IZhang. Cui & Chendfl997l) . Later in this 
paper we discuss the effects of the disc with the continuous torque 
through the inner radius. 

THCOMP is a thermal Comptonization model, using an 
approximate solution of the [Kompaneetsl ll956l) equation 
(Zdziarski, Johnson & Magdziarz 1996). It is parameterized by the 
asymptotic power-law photon index, T, electron temperature, T c 
and seed photon temperat ure, r see d- COMPPS is another Comp- 
tonization model IPoutanen & Svenssonlll996l) . It finds an exact 
numerical solution of the Comptonization problem explicitly con- 
sidering successive scattering orders. It also allows for various hot 
plasma geometries. For a slab geometry of the emission region, 
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Table 2. Summary of spectral models used in this paper. 



Model XS PEC components Description 



DTH 


WABS 


*(DISKBB+THCOMP) 


Multicolour disc and thermal Comptonization of the disc photons. 


DTF 


WABS 


*(DISKBB+THCOMP) 


The same as DTH, but the seed photon temperature is free and independent of 
the disc. 


DBTH 


WABS 


*(DISKBB+BBODYRAD+THCOMP) 


Multicolour disc, single-temperature blackbody and thermal Comptonization of 
the blackbody photons. 


DBTF 


WABS 


*(DISKBB+BBODYRAD+THCOMP) 


The same as DBTH, but the seed photon temperature is free and independent of 
the blackbody. 


DBPS 
DBPF 


WABS 
WABS 


*(DISKBB+BBODYRAD+COMPPS) 
*(DISKBB+BBODYRAD+COMPPS) 


The same as DBTH, with THCOMP replaced by COMPPS. 
The same as DBTF, with THCOMP replaced by COMPPS. 



this model is parameterized by the optical depth of the scatter- 
ing medium, r, electron temperature, T c , seed photon temperature, 
T scc d, and the inclination angle to the slab normal. The interstel- 
lar absorption was described by model WABS with the hydrogen 
column density TVh as a parameter. 

4.2 XMM-Newton 

iMiller et all 120031) successfully fitted EPIC-pn spectrum of 
XTE J 175 1-305 with a simple phenomenological model of a black- 
body and power law. We confirmed that this model gave a very 
good reduced x 2 /v = 2013/1859, indeed. There are, however, two 
fundamental issues here. Firstly, the most probable radiative pro- 
cess that gives rise to emission above ~1 keV is thermal Comp- 
tonization which can be described by a power law in a narrow en- 
ergy band only. At energies around seed photon temperature, the 
Comptonization spectrum has a low-energy cutoff. Secondly, the 
best-fitting spectral index of the power law (T = 1.44 ± 0.01) is in 
strong disagreement with RXTE data, where the spectrum is much 
softer, T ~ 1.8-1.9 (see below). 

Therefore, we tried a physically motivated model. Instead of 
the power law we used THCOMP, assuming that the seed photons 
were from the blackbody component. This model gave rather poor 
fit with X\ l v = 2754/1859 and broad residuals implying wrong 
shape of the continuum. In the second attempt we allowed the 
seed photons for Comptonization to be independent of the black- 
body component. This fit was much better { X 2 /v = 2044/1858), 
but the large apparent area of the blackbody, Abb ~ 1400 (D/8.5 
kpc) 2 km 2 , ruled out its origin from the neutron star surface. In- 
stead, it rather originated from the accretion disc, so we replaced 
the blackbody by the multicolour disc (model DTF in Table|2j- The 
results are shown in Table|5| The seed photons for Comptonization 
(fcTseed = 0.681q'qi keV) were much hotter than the disc photons 
(fcTdisc = 0.381q'q2 keV) and they might have ori ginated from 
the neutron star surface (see lGierfi hski & Done 2002; GDB02). We 
tried to model the surface emission by adding a single-temperature 
blackbody component with its temperature tied to the seed photons 
(see Fig.|2ji). Model DBTH improved the fit by Ax 2 = 28 over 
model DTF, with one degree of freedom less (Table[3J. The appar- 
ent area of this blackbody, Abb = 36^ (D/8.5 kpc) 2 km 2 , was 
consistent with a hotspot on the neutron star surface. 

There were no distinct features around 7 keV in the residuals, 
and we did not detect Compton reflection at a statistically signif- 



Table 3. Best-fitting parameters of the models applied to EPIC-pn data. 
Models are described in Tablel2l Inner disc radius i?; n and apparent area of 
the blackbody Abb were calculated assuming distance of 8.5 kpc. THCOMP 
temperature was fixed at 30 keV. 



Model DTF DBTH 



N n (10 22 cm" 2 ) 


1 04+ - 02 


1 ni+ 001 


fcT disc (keV) 




51+ 09 


i?i n V cos i (km) 


16 5 +25 


10.0+11 


fcTbt (keV) 




089lg;» 


Abb (km 2 ) 




3611 


r 


i.9o±r 01 


, 7? +0.07 
1 - /z -0.14 


kT sccd (keV) 


0.68«;° 2 


= fcTbb 


X 2 /" 


2039.5/1858 


2011.9/1857 



(a) 









^seed ^bb 



(b) 


*; 




Tbb 


T 

seed 


?bb 



Figure 2. Schematic geometries of the emission region, consisting of the 
accretion shock above, and the hotspot on the surface of the neutron star, 
(a) Geometry corresponds to a small temperature gradient in the hotspot, 
i.e. the case where seed photons have the same temperature as the pho- 
tons in the observed black body component (models DBTH, DBPS). (b) 
Geometry corresponds to a large temperature gradient where seed photons 
are hotter than the illuminated region of the neutron star surface around the 
shock, which gives rise to the observed black body emission (models DBTF, 
DBPF). The emitting region corresponds to a point where the accretion col- 
umn impacts on the surface (see fig. 12 in GDB02). 
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Figure 3. EPIC-pn/PCA cross-calibration. The figure shows the residuals 
(data-to-model ratio) of the EPIC-pn (triangles) and PCA (squares) spectra 
to a power-law model. The EPIC-pn spectrum was rebinned for presentation 
only. 



icant level in the EPIC-pn spectrum. The best-fitting amplitude of 
reflection was only Q/2tv = O.Oljlo'.oi giving fit improvement of 
A\ 2 = 1.9 with two degrees of freedom less. 



4.3 Simultaneous XMM-Newton and RXTE 

Encouraged by the results from XMM-Newton we extended the ob- 
served bandwidth by adding an RXTE spectrum. Observation 6 
from Table Q was simultaneous with EPIC-pn data. It overlapped 
with a small fraction of the EPIC-pn exposure only, but there was 
very little spectral variability throughout this observation, so spec- 
tra from both instruments should be the same. We also checked 
the cross-calibration between the two instruments to look for pos- 
sible response discrepancies. We fitted both spectra in the overlap- 
ping 3-10 keV band by a power law. The spectral indices were 
in good agreement, T = 1.660 ± 0.005 and 1.667 ± 0.015 for 
EPIC-pn and PCA, respectively. The residuals were also quite sim- 
ilar (see Fig. [5}, indicating nearly identical spectral shape. We 
decided that the broad-band analysis of the simultaneous XMM- 
Newton and RXTE was feasible. We fitted the EPIC-pn spectrum 
and PCA/HEXTE observation 6 together with a sequence of mod- 
els, which are summarized in Table [2] We allowed for the rela- 
tive normalization of these instruments to be free and used PCA 
normalization for flux calculations. The relative normalizations of 
EPIC-pn, HEXTE clusters and 1 with respect to PCA were 0.64, 
0.57 and 0.74, respectively, for our best fit. The fit results are shown 
in Tableland two of the fitted spectra in Fig. [4] 

First, we tested the possibility that the seed photons for Comp- 
tonization came from the accretion disc (model DTH). The model 
fitted the data fairly well (reduced x 2 l v = 1-15). It required how- 
ever a very small inner radius of the disc, Ri n V cos i = 3.6 ± 0.2 
(D/8.5 kpc) km (here i is the disc inclination angle), smaller than 
the expected neutron star radius. 

The model DTF, where the seed photons for Comptonization 
were independent of the disc temperature, gave a good fit to the 
EPIC-pn data in Section l4~2l We applied the same model to the 
broad-band spectrum (with free T c ). It improved the fit by A^ 2 = 



63.7 (at 1983 d.o.f.) with respect to DTH. The disc temperature 
was much lower and its inner disc radius much larger, making it 
consistent with an overall picture of a millisecond accreting pulsar 
(GDB02) or, more generally, an atoll source in the island state (e.g. 
lBarrel200ll) . 

Adding the blackbody source of the seed photons (presum- 
ably neutron star surface) as an explicit component (model DBTH) 
improved the fit further by Ax 2 = 53.8 (at 1982 d.o.f.). The ap- 
parent area of the blackbody component, Ahb — 49 ± 4 (D/8.5 
kpc) 2 km 2 , was consistent with a hotspot on the neutron star sur- 
face. However, this component did not produce enough photons for 
Comptonization. The Comptonization spectrum of T ~ 1.8 and 
kT c ~ 25 keV requires a similar luminosity in the seed and Comp- 
tonized photons, while the observed blackbody component was ~ 
30 times weaker than THCOMP. This means that majority of the 
seed photons were not visible, e.g. because they were covered by 
the Comptonizing accretion column (Fig.^). In fact, spectral fits 
with COMPPS (see below) showed that the actual area of the seed 
photons source greatly exceeded the expected area of the neutron 
star. 

Therefore, we considered one more model (DBTF), which in- 
cluded the disc emission, the blackbody and Comptonization of 
unseen seed photons of much higher temperature (see Fig. |2j)). 
This model produced a very good fit, improving previous results 
by Ax 2 = 27.3. It gave high seed photon temperature of fcT ace d = 
1.75±0.18keVandmuchsofterspectralindexofr = 1.94±0.05. 
As we show below in this section, the seed photons in this model 
were consistent with the hotspot on the neutron star. 

Next, we tested how sensitive the above results were to the 
choice of a particular model. We replaced THCOMP by another 
Comptonization code, COMPPS. It computes the exact numerical 
solution of the radiative transfer equation for Comptonization in a 
given geometry and allows finding the normalization (or apparent 
area) of the blackbody seed photons. We chose a geometry in which 
a slab of hot plasma is irradiated from below by the blackbody seed 
photons. The parameter r in the model is the vertical optical depth 
of the slab, so the actual line-of-sight optical depth depends on 
the inclination angle i. Unlike THCOMP, the unscattered seed pho- 
tons transmitted through the slab are already included in COMPPS 
model, so any additional blackbody component in the model cor- 
responds to photons not entering the hot plasma. This geometry is 
consistent with an accretion shock above the hotspot on the neutron 
star surface (GDB02; PG03). 

For this model we assumed that the inclination angle of the 
slab with respect to the observer was equal to the inclination of the 
system, i — 60°. This is, certainly, an approximation, as the ac- 
cretion shock rotates with the neutron star surface. However, as we 
show later in this paper (see Section l8~3l and Fig. 1151 the magnetic 
inclination 8 is small, < 10°, so the angle at which we see the shock 
doesn't vary much. 

In the model DBPS, corresponding to the geometry in 
Fig.|2ja), the temperature of the seed photons was tied to the black- 
body component. The apparent area of the source of the seed pho- 
tons was very large, A sccd = 790±l° 40 (D/8.5 kpc) 2 km 2 . This 
definitely could not be a hotspot on the surface, unless the distance 
to the source is less than ~3 kpc (but see Section l8^4l . 

A fit with seed photons independent of the blackbody tem- 
perature (model DBPF corresponding to a geometry in Fig. |2j>) 
gave better x 2 > but more importantly, yielded much smaller ap- 
parent area of the seed photons of 20!^ (D/8.5 kpc) 2 km 2 . On 
the other hand, the apparent area of the (independent) blackbody 
component was significantly larger, Abb ~ 90 (D/8.5 kpc) 2 km 2 . 
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Figure 4. Broad-band spectrum of XTE J1751— 305 from EPIC-pn, PCA and HEXTE. The models consists of the multicolour disc (green curve), single- 
temperature blackbody (blue curve), and thermal Comptonization of the blackbody photons (magenta curve). The dotted curves show unabsorbed spectral 
components. Panel (a) shows model DBTH, corresponding to the shock geometry shown in Fig.|5J. Panel (b) shows model DBTF, with geometry as in Fig.l2*b. 
See Table l2lf or model description and Tablel4lfor best fitting parameters. The lower panels shows the residuals of EPIC-pn (cyan), PCA (black) and HEXTE 
(blue). 



Table 4. Best-fitting parameters of several models applied to EPIC-pn/PCA/HEXTE simultaneous data. The models are summarized in Table l2l Inner disc 
radius i?; n and apparent area of the blackbody A bb are calculated assuming distance of 8.5 kpc. The spectrum and components of the model DBTH are shown 
in Fig.|] 
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The great advantage of COMPPS model is that it computes 
angle-dependent spectrum. For an intermediate optical depth of the 
slab of r ~ 1.5-2.0 the overall spectral shape of the Comptonized 
continuum depends only weakly on the observer's angle with re- 
spect to the slab. The main effect is the amount of unscattered seed 
photons seen at different angles, as it is very sensitive to the actual 
optical depth along the line of sight, r/ cos i. So far, we have as- 
sumed quite large inclination angle of 60° . To check the angular de- 
pendence we repeated the fits with models DBPS and DBPF at 30°. 
The DBPS fit was dramatically worse by Ax 2 = 75.4. Smaller an- 
gle gave larger, not consistent with the observed spectrum, amount 
of unscattered seed photons. The DBPF fit, with more freedom of 
shaping the low-energy part of the spectrum gave the fits worse by 
Ax 2 = 7. 1 and optical depth of r = 2.54 ± 0.15, larger than at 60° 



by expected factor of cos 30°/ cos 60°. The actual line-of-sight op- 
tical depth to the seed photons was about 2.9, so the observed un- 
scattered fraction of the seed photons was about 5 per cent. This 
explains why THCOMP, with no unscattered blackbody photons in- 
cluded in the model, could fit the data well. When inclination was 
allowed to be free, the best-fitting values were i = 79+35 and 
65tjg deg, for DBPS and DBPF, respectively. 

No reflection was significantly detected in any of these mod- 
els. We found only upper limits on Q/2n of 0.06 and 0.41, from 
DBTF and DBPF, respectively. The first model incorporated reflec- 
tion with a self-consistent iron line IZvcki. Done & S mith 1998), 
while the second one included an independent DISKLINE model, 
with one more free parameter and a less strict constrain. 

The total unabsorbed bolometric flux from model DBTF was 
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Figure 5. Fitting results of the blackbody plus thermal Comptonization 
model (DBTF from Tablel2l without disc component) to the outburst RXTE 
spectra. The electron temperature was fixed at 42 keV, the absorption col- 
umn at 1.0 X 10 22 cm -2 . T is the asymptotic photon spectral index, 
fcT seet j (keV) is the seed photon temperature, -/Vh ar d (10 -9 erg s — 1 cm -2 ) 
is the normalization of the Comptonized component at 10 keV, fcT so f t 
(keV) and A so f t (km 2 ) are the temperature and apparent area of the soft 
(blackbody) component, respectively. The fraction of Eddington luminos- 
ity, ibol/^Eddi and ^aoft are calculated for a distance of 8.5 kpc. 



Fboi = 2.9 x 10 erg s~ cm - , which corresponds to luminosity 
Lboi = 2.5 x 10 37 erg s" 1 or 0.12Z/Edd for a 1.4 M© neutron star 
at a distance of 8.5 kpc. The contribution to the total luminosity 
from the disc and blackbody components was 9 and 13 per cent, 
respectively. 



4.4 RXTE outburst 

Finally, we analysed the RXTE spectra from Table Q covering the 
whole outburst (observations 2-21). For spectral fitting we used 
some of the models from Table [2] however without the disc com- 
ponent, which could not be constrained by the PCA data and had 
little effect above 3 keV. We also fixed the absorption column at 
Nh = 1.0 x 10 22 cm -2 . 

Model DBTH gave a good total (i.e. summed over all 20 ob- 
servations) x 2 /v = 2153.7/2440. The spectral parameters remained 
roughly constant throughout the outburst. In this model the soft 
blackbody component was very weak, its normalization (or appar- 
ent area) consistent with zero in about half of the fits and its contri- 
bution to the bolometric luminosity negligible. On the other hand 
its temperature was well constrained, as it was equal to the tem- 
perature of the seed photons in the Comptonized component and 
constrained by the low-energy cutoff in the spectrum. 

When we made seed photon temperature a free parameter 
(model DBTF), some of the fits become unstable, in particular the 



electron temperature was not constrained. We fixed kT c = 42 keV 
(the value taken from the outburst average, see below). After this 
modification we obtain total x 2 l v = 2106.4/2240. The results are 
presented in Fig. [5] The blackbody and seed photon temperatures 
were slightly higher and the spectral index slightly softer in the 
beginning of the outburst, but apart from that the spectral shape re- 
mained amazingly constant throughout the outburst. The ratio of 
unabsorbed fluxes from Comptonization and blackbody was ~ 6. 

We also fitted the average spectrum of the entire outburst 
(co-added observations 2-21). Model DBTF gave a very good fit, 
X 1 l v = 101.4/121. The blackbody temperature was fcTbb = 
0.83 ± 0.06, seed photon temperature fcT sccd = l-46± ;i| keV, 
the spectral index T — 1.93 ± 0.03 and electron temperature 
kT c = 42I7 keV. These parameters were consistent with those ob- 
tained from the combined XMM and RXTE fits. Compton reflection 
was not significantly detected again: the was only an upper limit on 
Q/2n < 0.15, with fit improvement of Ax 2 = 2.0 and two degrees 
of freedom less. This average outburst spectrum formed a base for 
analysis of the phase-resolved spectra in Section[6| 



4.5 Summary of spectral results 

In the above section we have fitted the X-ray spectrum of 
XTE J175 1-305 by models consisting of the emission from ac- 
cretion disc, blackbody hotspot and optically thin Comptoniza- 
tion. We have considered two types of models, corresponding to 
hotspot/shock geometries sketched in Fig. [2] In the first geometry 
the seed photons for Comptonization originated from the hotspot 
(DBTH and DBPS, Fig. [4^). In the second geometry the seed 
photons were considerably hotter than the blackbody (DBTF and 
DBPF, Fig.|4j>). The second geometry was preferred both from sta- 
tistical point of view (it gave better fit) and from the physical con- 
straints, as the first geometry required the area of the seed photons 
source larger than the neutron star surface. In the geometry of the 
short shock above the surface, GDB02 ruled out the origin of the 
seed photons from the accretion disc. 



5 PULSE PROFILES 

Fig.|6|shows the pulse profile of XTE J 175 1-305 accumulated over 
the entire outburst. It was not perfectly symmetric, a fit by a single 
sine function gave rather poor x 2 l v — 28.3/13. With the second 
harmonic added (its phase independent of the first one), the fit was 
much improved, x 2 l v — 17.3/11. The rms amplitudes, calculated 
from these fits, were 3.28±0.03 and 0. 11±0.03 per cent for the 
first and second harmonic, respectively, which correspond to the 
peak-to-peak amplitudes (F max -F min )/(F max + F min ) of 4.6 and 
0.15 per cent. The lower panel of Fig.|6|shows the residuals of the 
single-harmonic fit with the second harmonic function overplotted. 

We investigated the energy dependence of the pulse profile. 
Fig. Q shows the outburst-averaged pulse profile in three energy 
bands. The pulse amplitude decreased and its profile shifted to- 
wards earlier phase with increasing energy. This effect can be 
clearly seen in Fig. [8] where we calculated the amplitude and phase 
lag in each of the PCA energy channels by fitting their pulse pro- 
files by a sine function. Fig.|8ji shows the dependence of the pulse 
rms amplitude on energy, which decreased from around 3.3 per cent 
at 2-10 keV to less than 2.5 per cent at 20 keV. The absolute value 
of the time lag with respect to the energy channel 2-3.3 keV in- 
creased with energy up to about 120 /is at 10 keV, above which it 
seemingly saturated (Fig.lSJj). The negative time lag means that the 
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Figure 6. Pulse profile of XTE J175 1-305. The upper panel shows the 
folded 2-60 keV PCA light curve accumulated over the entire outburst. The 
eiTors are of order of symbol size. The curve represents the model consist- 
ing of two harmonics. The lower panel shows the residuals of data minus 
first harmonic and the second harmonic model curve. 
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Figure 7. Pulse profiles (outburst averaged) in three energy bands. 



hard photons arrived earlier than the soft photons. This behaviour 
was very similar to other millisecon d pulsars, SAX J1808.4— 
3658 ICui. Morgan & Titarchukl fi998: GDB02), XTE J0929-314 
iGallowav et alJ l2002h . and IGR J00291+5934 iGallowav et ail 
l2005t) . 



6 PHASE-RESOLVED ENERGY SPECTRA 

The photon statistics in the phase-resolved spectra was very lim- 
ited, as each of the phase bins contained only a sixteenth part of 
the total counts. This made spectral fitting of individual observa- 
tions difficult. However, as we have shown in the previous section, 
the spectral shape remained almost constant throughout the out- 
burst, so analysis of the outburst-averaged data was feasible. Thus, 
we fitted the average phase-resolved spectrum in each of the phase 
bins by the models developed earlier in this paper. 

The best physically motivated models for the broad-band data 
were DBTF and DBPF, from which we chose DBTF as it is much 
faster to compute. We simplified this model by removing the disc 
component, which is not constrained by the PCA-only data. We 
fixed N H = 1.0 x 10 22 cm -2 and kT c = 42 keV. This, however, 
proved to be problematic. Though the total \ 2 l v (i- e - summed over 
all 16 phase bins) was very good (407/672), there was a strong 



Figure 8. (a) The rms amplitude of the pulse profile as a function of energy, 
(b) Time lag of the pulse profile versus energy, with respect to the first 
energy channel, 2-3.3 keV. 



correlation between the spectral index and seed photon tempera- 
ture in the limited PCA bandpass. When both parameters were al- 
lowed to be free, the fits tended to give very low T see d, hard F, 
and much higher T so f t than in the broad-band fits. These results 
were inconsistent not only with the EPIC-pn data, but also with 
the PCA/HEXTE phase- and outburst-averaged spectrum. There- 
fore, we forced softer spectral index, and fixed T — 1.93, the value 
obtained from the PCA/HEXTE outburst-average fit (Section l4.4t . 
This gave worse overall x / v - 480/688, but the fitting parameters 
were consistent with broad-band fits. 

We followed GDB02 in the fitting procedure to eliminate spec- 
tral parameters that did not pulsate significantly, i.e. for each the 
intrinsic variance [eq. (1) in GDB02] was zero. In each step of the 
fitting procedure we fixed one of the non-pulsating parameters at 
its mean value from the previous fit, until we ended up with just 
two pulsating parameters: normalization of the blackbody (its ap- 
parent area) and normalization of THCOMP. The total x 2 l v was 
501/720. The result is presented in Fig.|5| The normalization of the 
Comptonized component can be well fitted by the profile 



A r hard(0) = 0.41 + 0.017 sin[27r(0 - 0.21)] 10 



-9 —1 —2 

erg s cm , 



where <j> is the phase. The blackbody apparent area can be described 
by 

A so{t ((p) = 101 + 6.4sin[27r(0 - 0.30)] km 2 . 

The intrinsic variance of iVhard and A BO f t was 3.1 ± 0.1 and 
4.4 ± 0.4 per cent, respectively. The second harmonic was not sig- 
nificantly present in any of these profiles, with upper limits of 8 and 
12 per cent of the first harmonic, in /Vhard an d -Asoft. respectively. 

Next, we checked if the spectral variability with phase can be 
reproduced not by two independently pulsating components, but by 
change in the spectral index of Comptonization. We modified the 
above spectral fits, starting with the model where F was free but 
the seed photon temperature was fixed kT scc d ~ 1-46 keV (Section 
14.41 . In the last-but-one step three parameters were allowed to vary: 
^4soft, A?hard and r. Both normalizations pulsated with non-zero 
intrinsic rms, while the intrinsic rms of T was zero. When we fixed 
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Figure 9. Pulsation in the model consisting of the blackbody and Comptonization of the seed photons independent of this blackbody. Only normalizations of 
both components were free during the fits. This model with total x 2 / u = 501/688 is preferred over the pulsating index model from Fig. 1101 




Figure 10. Pulsation in the model consisting of the blackbody and Comptonization of the seed photons independent of this blackbody. Only normalization of 
the hard component and its spectral index were free during the fits. This model gave worse fit, x 2 l v — 592/720, than the two pulsating components model 
from Fig.l9l 



^4 so ft and allowed only iVhard and T to vary, we obtained the overall 
fit with x 2 /v = 592/720, much worse than in the previous case 
(see Fig. HOt . Clearly, the model with two components pulsating 
independently was preferred over the pulsating spectral index. 

In reality, these two approaches may not be sepa- 
rated. In Comptonization models, the angular distribution of 
the escaping radiation depends on the scattering order (see 
e.g. ISunvaev & Titarchuklll985t Iviironen & Poutanenl l2004 and 
Fig. I16> . Thus the spectral index of Comptonized radiation is a 
weak function of the viewing angle, therefore, a slab of Comp- 
tonizing plasma on the surface of the spinning star can produce a 
spectrum with pulsating spectral index (as in the second approach). 
At low energies, one expects a contribution from the seed, black- 
body photons which, being a strong function of the inclination 
exp(— t I cos i), can give rise to the pulsations of relative normal- 
izations of soft and hard spectral component (as in the first ap- 
proach). To check this, we have fitted the phase-resolved spectra 
with COMPPS model, which produces angle-dependent spectrum, 
where two parameters, inclination and normalization, were allowed 
to vary. The fits occurred to be rather unstable, producing multiple 
minima in % 2 f° r small and large inclinations. The high- angle so- 
lution yields the mean slab inclination of about 68° and amplitude 



of only about 2° . Such a small variation of i would be possible if 
the shock (and magnetic pole) were almost perfectly aligned with 
the rotational pole of the star. 



6.1 Summary of phase-resolved results 

In the section above we have applied two distinctly different ap- 
proaches to the phase-resolved spectra. In the first approach there 
were two components of fixed spectral shape (blackbody and 
Comptonization), varying independently with phase. In the second 
approach the spectral slope of the hard component (Comptoniza- 
tion) varied as a function of phase. The crucial common feature 
of these two models was the lag in pulsation between softer and 
harder energies. This manifested itself as energy-dependent time 
lags (Fig. [SJ)) where the hard photons arrive earlier than the soft 
photons. The model with two independently varying components 
was preferred over the varying index. In reality, both models can 
operate together. 
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Figure 11. Power density spectrum of observation 2 from TablelTlfitted by 
a multi-Lorentzian model. 



7 POWER DENSITY SPECTRA 

Multiple Lorentzian profiles can give a good description of the 
power density spectra (PDS) of black hole (e.g. Nowak 2000) and 
neutron star binaries (e.g. van Straaten et al. 2002). We followed 
this approach and fitted the PDS of all observations from Table Q 
by a model consisting of up to three Lorentzians. To avoid the un- 
certainties of the high-frequency part of the spectrum we limited 
the fit to the frequency band of 0.01-100 Hz. An example of the 
PDS with the best-fitting model is presented in Fig. II II 

The low- and high-frequency Lorentzians, describing the 
broad noise components were zero frequency-centred. The mid- 
frequency Lorentzian function had its centre frequency free. The 
characteristic frequency (width) of the low-frequency Lorentzian, 
i^lf, can be attributed to the low-frequency break in the PDS. The 
centre frequency of the mid-frequency Lorentzian, vqpo, corre- 
sponds to the quasi-periodic oscillation (QPO) at 2-6 Hz. Evolu- 
tion of these two frequencies throughout the outburst is shown in 

Fig.m 



8 DISCUSSION 

8.1 Comparison with SAX J1808.4-3658 

GDB02 built a physically-motivated model of another millisec- 
ond pulsar, SAX J1808.4-3658, analysing its phase-averaged and 
phase-resolved X-ray spectra. According to this model the accre- 
tion disc is disrupted at some distance from the neutron star by its 
magnetic field, and collimated towards the poles to create an accre- 
tion column impacting on to the surface. The material in the col- 
umn is heated by a shock to temperatures of ~40 keV. The hotspot 
on the surface with temperature < 1 keV provides seed photons 
for Comptonization in the hot plasma. Due to a different character 
of emission from the optically thick spot and optically thin shock 
there is a phase shift between the soft and hard photons. PG03 ex- 
tended this model to incorporate relativistic effects and different 
angular distribution of emitted radiation from the hotspot and the 
shock. They showed that the emission patterns expected from the 
optically thick blackbody and optically thin Comptonization can 
explain the observed phase-resolved spectra indeed. 



Figure 12. Evolution of characteristic frequencies in the PDS through- 
out the outburst. The open squares correspond to the width of the low- 
frequency, zero-centred Lorentzian, z^lf- The filled squares describe the 
QPO frequency, fQpr> 



The X-ray properties of XTE J175 1-305 are remarkably sim- 
ilar to SAX J1808.4-3658. Th e outburst light curv e has an alike 
exponentially decayed profile I Gilfan ov et alJll998h with a break 
or re-flare iShahbaz. Charles & Kina 1998i) about 13 days into 
the outburst. The 3-200 keV spectral sh ape, entirely typical of an 
atoll source in the hard/i sland state (e.g. lBarret & Vedrenndfl994l 
iGierlinski & Dondl2002h . is very similar to SAX J1808.4-3658, 
and l ikewise almost constant during the outburst iGilfanov et alJ 
1998; GDB02). Both sources showed coherent pulsations at 
around 400 Hz, with similar rm s amplitude of a few per cent 
fWiinands & van d er Kli j I l998al) and similar energy-dependent 
soft time lags <Cui et alJll99 Si. The power density spectra below 
100 Hz have similar shape and frequencies of the QPO and break 
in the s pectra correlated in a similar way l Wiinands & van der Klis 
HHH). 

The peak outburst luminosity of XTE J175 1-305 was about 
an order of magnitude greater than in SAX J1808.4-3658, but the 
distance to the former source is uncertain (see discussion in Section 
lOt . 

The main difference in spectra between SAX J1808.4-3658 
and XTE J175 1-305 was the weaker blackbody component in the 
latter source. This can be easily explained by larger optical depth of 
the hot plasma, if the blackbody hotspot is covered by the Comp- 
tonizing shock. As GDB02 used models in different plasma geom- 
etry than in this paper, it was difficult to compare them directly. 
Therefore, we have fitted the outburst average PCA/HEXTE spec- 
trum of each source with the same model consisting of a black- 
body and its Comptonization (COMPPS) in slab geometry, inclined 
at 60° (corresponding to DBPS model without the disc). We found 
kT c = 60±g keV and r = 0.88l ;^ for SAX J1808.4-3658 and 
kT c = 33±2 keV and r = 1.7 ± 0.1 for XTE J1751-305. Tak- 
ing into account the inclination of the slab, the line-of-sight op- 
tical depth was 1.8 and 3.4, and the amount of unscattered seed 
photons reaching the observer was 17 and 3 per cent, respectively. 
One can note here that the product rT c is almost the same in the 
two sources. This is natural in the two-phase models (cool neutron 
star plus a hot shock above), where the temperature depends on the 
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optical depth, because of the ener gy balan ce constraints (see e.g. 
Haardt & Maraschj|l993l;IStern et alJ| l995tlPoutanen & Svenssonl 
ll996tlMalzac. Beloborodov & Poutanenl200lh " 

The pulsation of the soft component in XTE J 175 1-305 was 
significantly smaller (intrinsic variance of 4.4 per cent) than in 
SAX J1808.4-3658 (18.1 per cent), though we must stress that dif- 
ferent models were used to estimate these numbers. The pulse pro- 
file of the hard component was less skewed than in SAX J1808.4— 
3658. The upper limit on the second harmonic was 8 per cent of 
the first harmonic, while SAX J1808.4-3658 required the second 
harmonic at 18 per cent level. The relative weakness of the sec- 
ond harmonic in XTE J1751-305 is expected, since the ratio of 
the amplitudes of the second to the first harmonic is proportional 
to the product sin i sin 6 iPoutanenlEool . where i is the inclina- 
tion of the rotational axis and 8 is the co-latitude of the hotspot. At 
the same time, the variability amplitude is also proportional to the 
same product (see eq. (8) below). Thus the smaller is the observed 
variability amplitude, the closer the profile is to a sinusoid. 

Phase lags between the soft and hard photons were also 
small er in XTE J175 1—305 (see Fig. [SJs in this paper and fig. 4 in 
ICuietalJl9 98). In both sources the time lags are negative, i.e. hard 
X-rays arrive before soft X-rays. The energy dependence was also 
similar. The absolute value of phase lags increased with energy up 
to ~10 keV, above which they saturate at different levels: ~ 200 /is 
(8 per cent of the pulse period) for SAX J1808.4-3658 and ~ 100 
fis (4 per cent) for XTE J175 1-305. Smaller values of the phase 
lags can result from a weaker soft component - the black body 
has a very different radiation pattern from the Comptonize d emis- 
sion r esulting in a different pulse profile (see lViironen & Poutanenl 
2004). 

There was no Compton reflection significantly detected in 
any of the models used in this paper, while GDB02 reported 
weak (Q/2n ~ 0.1) but statistically significant reflection in 
SAX J1808.4-3658. On the other hand, the upper limit on the re- 
flection amplitude from RXTE data only (as used by GDB02) was 
0.15, making it consistent with SAX J1808.4-3658. 



8.2 Spectral components 

8.2.1 Origin of the soft components 

The best-fitting model to the broadband data required three com- 
ponents (Fig.0J: two soft components ('blackbodies') and the hard 
one ('power law'). We interpreted them as emission from the ac- 
cretion disc (low-temperature soft component), the hotspot on the 
neutron star surface (high-temperature soft component) and Comp- 
tonization in the shock above the stellar surface (the hard compo- 
nent). 

The simplest argument for this interpretation of soft compo- 
nents came from their apparent area. We can easily compare them 
directly by replacing the disc component in DBTF model with an- 
other single-temperature blackbody, so the model would consist of 
two blackbodies and Comptonization. We fitted this model to the 
broadband spectrum and obtained a good fit. As a result, the appar- 
ent area of the low-temperature blackbody was ~10 times larger 
than the hotter one. 

Furthermore, we would not expect any pulsation from the disc, 
as it rotates independently of the neutron star. Since the reflection is 
low, irradiation of the disc by the spinning spot/shock is negligible. 
The soft component in the PCA band (the hotter soft component 
from Fig.[4j showed clear pulsation with about 4 per cent rms, so 
it did not originate from the disc. We can also rephrase this argu- 



ment in a model-independent way. If the disc emission had strongly 
contributed to the PCA spectrum, we would see drop in the pulse 
rms at soft X-rays, but it was not the case (Fig. [8^). Therefore, the 
disc emission must be weak above 3 keV, and the observed hotter 
blackbody rather comes from the neutron star. 



8.2.2 Origin of the hard X-ray emission 

Since the hard X-ray emission is pulsed, a fraction of it must origi- 
nate from regions confined by the magnetic field. The most obvious 
source of hard X-rays is the place where material collimated by the 
magnetic fields impacts on to the surface iBasko & Sunvaevll976l) . 
However, the pulsed emission can also be released in the m agne- 
tosphere further away from the surface (e.g. lAlv & Kuiipersl 1990). 
In addition, an unknown fraction of th e X-rays can come from t he 
corona above the accretion disc (Galeev, Rosner & Vaiana 1979). 

The X-ray spectrum of XTE J175 1-305 was fairly constant 
as a function of pulse phase and our best-fitting model consisted 
of two components of constant spectral shape pulsating indepen- 
dently. This favours geometries where all of the emission comes 
from the same source, as opposed to the sum of separate constant 
and pulsed components. The idea of the bulk of emission origi- 
nating from polar caps is also supported by the observed broad- 
ening of the pulse peak in the PDS in SAX J1808.4-3658 due 
to mod ulation of the aper iodic variability by the harmonic spin 
period iMenna et alll2003t) . We cannot yet confirm this result for 
XTE J175 1-305, but the analogy between the two sources makes a 
strong argument for the emission predominantly coming from the 
short, heated shock in the polar region. 

The constancy of the spectral slope during the whole outburst 
can be used as an argument that the emission region geometry does 
not vary much with the accretion rate. If the energy dissipation 
takes place in a hot shock, while the cooling of the electrons is 
determined by the reprocessing of the hard X-ray radiation at the 
neutron star surface (two-phase model, Haardt & Maraschi|^9S 
IStern et all 1 9951: iPoutanen & Svenssonl 1996tlMalzac et all200ll) 
the spectral slope is determined by the energy balance in the hot 
phase and, therefore, by the geometry. 

Our spectral models in Section 14751 were consistent with the 
slab geometry. For the two models with COMPPS we assumed a 
slab of hot plasma irradiated from below by a blackbody. The seed 
photon area can be determined from the Comptonization theory be- 
cause the amplitude of the Comptonized component depends lin- 
early on the seed photon flux. In model DBPS (see Table |4j the 
seed photons came from the blackbody component of the spectrum, 
with temperature of ~ 0.8 keV. This would correspond to the emis- 
sion region geometry depicted in Fig. |2a). However, the area of 
the seed photons in this model, i.e. the spot under the slab, was 
very large, j4 scc d ~ 800 km 2 , for a distance of 8.5 kpc (see discus- 
sion of distance in Section[8~4). In the best-fitting model DBPF the 
seed photons were independent of the blackbody component. This 
could correspond to a spot with a hotter centre and cooler edges 
(Fig.|2p). The central part of the spot with temperature of ~ 2 keV 
was covered by the Comptonizing shock and was not visible di- 
rectly. This constituted the main source of the seed photons. The 
cooler (~ 1 keV) outer part of the spot could be seen directly. The 
apparent area of the inner and outer spot were A ~ 20 and ~ 90 
km 2 , respectively, assuming distance of 8.5 kpc. 
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Figure 13. Trajectory of a photon in the gravitational field of a neutron star. 
Photons emitted at angle a > a c relative to the stellar normal cross the 
disc plane. 



8.3 Geometry 

8.3.1 Inner radius of the accretion disc 

The spectral models from Section l431 predicted the disc to be trun- 
cated fairly close to the neutron star, at J?i n \/cos i > 10 km. Our 
preferred model DBPF yielded J?j n w 14 km for the inclination of 
60°. The inner disc radii from Table |4] were calculated assuming 
the torque-free inner boundary condition, so the disc temperature 
was: 



T e ff oc R- 




(1) 



The optically thick disc might be truncated into an inner 
hot flow, as it has been sugg e sted for black holes (e.g. 
I Shapiro. Lightman & Eardlevlll97a lEsin. McClintock & Naravanl 
1997), or disrupted by the star's magnetic field. In either case it 
is possible that t he torque is retained through the transition region 
(e.g. Agol & Krolik 2000). Dropping the inner boundary condition 
term from Eq. Q would have increased the derived Ri n by fac- 
tor ~ 2.7 ( Gierl mski et alJl999T) . in which case DBPF would have 
given Ri n « 40 km. 

Another important factor is the uncertainty in the distance to 
XTE J175 1-305. The above values of R[ n were calculated for the 
distance of D = 8.5 kpc. If the source were at 3 kpc (see Section 
18.41 . we would have Ri n ~ 5 and 13 km from DBPF model with 
torque-free and torqued boundary condition, respectively. 

We can estimate the inner disc radius independently from the 
observed amplitude of reflection. Let us assume that the emission 
originates from a small spot close to the rotational pole of the neu- 
tron star. (These assumptions give a lower limit on reflection am- 
plitude.) Because of the light bending, a photon emitted at an an- 
gle a relative to the rotation axis, travels at the infinity along the 
trajectory which makes an angle ip with the axis (see Fig. 1131 - 
The relation between a and tj) i s given by an elliptical integral 
( Pech enick. Ftaclas & Cohenll983l) . At small a, photons propagate 
in the upper hemisphere never crossing the disc plane. A photon 
emitted at a critical angle a c flies parallel to the disc surface. A 



simple estimate for this angle can be obtained using Belo borodovl 
(2002, hereafter B02) approximation for light bending: 

cos a = R s /R, + (1 - Rs/R») cosip, (2) 

where R, is the stellar radius and R B = 2GM/c 2 is the 
Schwarzschild radius for a star of mass M. Substituting ip = tt/2, 
we immediately get cosa c = R s /R,. Photons, emitted at an an- 
gle a satisfying condition < cos a < Rs/R,, will cross the disc 
plane at a radius given by an approximate expression for the photon 
trajectory (B02) 



Rdis 



R* 



(1 - cos?//) 2 
4(1 + cos?/;') 2 ' sin 2 ?/;' 



+ 



1/2 



- cos 



2(1 + cos?/>') 



(3) 



where the impact parameter b — R, sina/^/l — R s /R,, ?/>' = 
ip — n/2, and if> is related to a via equation J2j. The smallest ra- 
dius where photon trajectories can cross the disc is obtained by 
substituting a — n/2, sin?/)' = — 1/(R,/R S — 1) and b — 
R,/\/l - Rs/R, to equation (5}. 

We now can compute the reflection amplitude as a function 
of the inner disc radius. Let us assume for simplicity that radiation 
from the spot has a constant specific intensity 1(a) — Jo = const 
(as for the black body emission). We can compute the luminosity 
(from the unit area) escaping to the infinity as 



Lout = 2ty Io cos a sin a da = 
Jo 



f. (4) 



If the inner disc radius is larger than J?di S c(7r/2), we reverse rela- 
tion |3} to find the emission angle a^ n at which a photon should be 
emitted to cross the disc at radius R[ n . The reflection luminosity is 
then 



2tt 



Jo cos a sin a da — 



Jo 



(5) 



If Rin < J?diac(7i"/2), the reflection luminosity is maximal 

Jrcfl,max = (Rs/R*) 2 Io/2. 

In order to estimate the reflection amplitude, we need to spec- 
ify the angular dependence of the observed escaping and reflected 
luminosities. The gravitational redshift affects both luminosities 
the same way and thus can be neglected. The escaping luminos- 
ity depends on the inclination angle as (see eq. J2]) L out (i) oc 
R s /R, + (1 — Rs/R*) cosi, where normalization factor is found 
from the condition L ut = / J/ ou t(i)d(cos i). We can assume that 
the reflected luminosity is proportional to the disc area projected on 
the sky: L Ten (i) = L re fl2cosi. Thus the reflection amplitude is: 



2^ 



Li, 



(1 + Rs/R*) cosi 



(6) 



Lout Rs/R* + (1 - Rs/R*)cosi ' 

It is plotted in Fig. ll4l for different stellar compactnesses and incli- 
nation angles. Because the observed upper limit on the reflection 
amplitude is Q./2-r < 0.07 (see Section l4~2l . the inner disc can- 
not be closer than about 90 km for inclination smaller than 60°. At 
i — 80°, constraints become much weaker and depend on stellar 
compactness. If stellar radius J?* < 2.5J? a = 10.5 km (neutron star 
mass of 1.4M0 was assumed), then Ri n > 40 km, but if J?* = 3J? S 
reflection is smaller (since bending is weaker) and any inner radius 
would satisfy these constraints. 

Menn aet alj J5)03) estimated the magnetospheric 
radius for SAX J1808.4-3658 to be ~17 km and ar- 
gued that the accretion disc is di srupte d at this radius. 
Ivan Straaten. van der Kris & Wiinandsl J2004I) analysed QPO 
frequency correlations in SAX J1808.4-3658 and found that 
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Figure 14. Reflection amplitude as a function of the inner disc radius. The 
source of photons is a small spot at the rotational pole of the neutron star. 
The intrinsic radiation pattern is assumed to be 1(h) = const (as for the 
black body emission). Dotted curves correspond to R* = 2R B , solid curves 
to R* = 2.5-Rs and dashed curves to R t = 3R B . The upper curves are for 
the inclination i = 30°, middle ones for i = 60°, and the lower curves are 
for i = 80°. Neutron star mass of 1.4 Mq is assumed. 



the upper kilohertz QPO has a frequency by factor ~1.5 lower 
when compared to other atoll sources, while low-frequency QPOs 
were comparable. They suggested that stronger magnetic field 
of SAX J1808.4-3658 moved the inner edge of the disc further 
away, and decreased the Keplerian frequency, responsible for the 
upper kHz QPO, while low-frequency QPOs were formed further 
out, and were not a ffected. There was no s imilar frequency shift 
in XTE J 175 1-305 Ivan Straaten et alj|2004) . suggesting weaker 
magnetic field and smaller disruption radius than aforementioned 
17 km. 

Our results about the inner disc radius are not conclusive. The 
direct disc emission observed in the broad-band spectrum corre- 
sponded to the disc truncated close to the neutron star. This result 
depended on the assumed disc model and a radius of ~ 40 km was 
possible for the disc without the torque-free inner boundary condi- 
tion. Similar or larger radius was estimated from the rather small 
amount of Compton reflection from the disc. If the inner disc ra- 
dius had been larger than the magnetospheric radius, it would have 
been truncated into an optically thin hot inner flow, similarly to the 
low/hard state of black holes l Shapiro et al. 1976; Esin et al. 1997). 
We would like to present an argument in favour of this scenario. 

The low-frequency QPO and a break in the power spectrum 
are likely to be related to the truncation radius of the disc. The 
shape of the low-frequency power spectrum and frequency cor- 
relati ons are very similar in b lack holes and neutron stars (e.g. 
Wiinands & van der Klisll999T) . so their origin should be indepen- 
dent of the neutron star surface and magnetic field. The lack of 
correlation between the break/QPO frequency and the luminos- 
ity (see Figs.fT1and ll2> supports this idea. If the disc were trun- 
cated at the magnetospheric radius indeed, we would expect its 
inner radius Ri n = i? magn oc L~ 2 ^ 7 , where L is the luminos- 



ity (e.g. lHavakaw3ll985l) . Let us assume that the low-frequency 
QPO is connected to the precession time-scale of a ver tical pertur- 
bation in the disc iStella & Vietrilll998tlPsaltis & Normarfe oOO). 
This yields vqt>o oc R^^ 2 , so we should expect z^qpo oc L 3 ^ 7 . 
Any other link between the QPO frequency and the truncation ra- 
dius would have resulted in the correlation between vqpo and L. 
This has not been observed. Therefore the optically thick disc was 
probably truncated above the magnetospheric radius, by a mech- 
anism independent of the magnetic field. This could be the same 
mechanism truncating the disc of the low/hard state in black hole 
binaries, e.g. evaporation (see Rozanska & Czernv 2000). The low- 
frequency part of the PDS is formed in this region. Below the trun- 
cation radius the disc is replaced by a hot, optically thin inner flow, 
which is eventually disrupted by the magnetic field at or near the 
magnetospheric radius. 



8.3.2 Spot size 

Let us estimate the physical size of the emission region using the 
model DBPF (hot shocked region surrounded by a cooler black 
body edges, see Fig. |2j>) which seems physically more realistic. 
The inner spot area of A ~ 20 km 2 corresponds to the 'observed' 
at the infinity radius of R^ = 2.5 km. These estimations, how- 
ever, do not take colour hardening and gravitational redshift and 
light bending into account. 

The observed blackbody (i.e. colour) temperature is differ- 
ent from the effective temperature by a colour hardening fac- 
tor, /coi = Tcoi/Tcff- For the hydrogen and helium atmo- 
spheres of weakly magnetized neutron stars, the exact value of 
the colour correction is not very well known, but it has been 
estimated to be / co i ^1. 3 jLewiikvan Paradiis & Tamml]P393j; 
IZavlin. Pavlov & Shibanolll996tlMadei. Joss & Rozariskdhool) . 
This would result in a blueshift of the blackbody spectrum, and the 
effective radius of the blackbody would be R B pat,eS = /coi^co- 
Spectral hardening results from the fact that the absorption opac- 
ity rapidly decreases with increasing photon energy and one sees 
deeper and hotter layers at high en ergies. Situation in acc reting 
pulsars can be quite opposite (see e.g. Deu fel. Dullemond & Spruit! 
2001): the hot layers are at the top and low-energy photons having 
larger absorption cross-section come from the hotter outer layers. It 
is possible that the resulting colour correction could be negligible, 

/col ~ 1. 

Implementing gravitational corrections requires quite compli- 
cated numerical treatment. However, in the case when the (small) 
spot is always visible (which is the case for XTE J1751-305, be- 
cause we observe almost sinusoidal oscillations) one can obtain a 
simple relation between the observed size at infinity and the phys- 
ical size at the neutron star surface i? spo t = -RooQ 1 ^ 2 (B02; 
PG03), where Q = R s /R, + (1 - R s /R*)cosi cos (9, and 6 is 
the colatitude of the spot centre (magnetic inclination). The small- 
est possible radius is obtained for i = 6 = 0, then Q = 1 and 
-Rspot = Roo- This corresponds to the minimum angular size of 
the spot of p = 12° (14°) for R, = 3R S (2.5i? s ). A reasonable 
upper limit can be obtained taking i = 90° and 8 = (or i = 
and 6 = 90°), now Q = Rs/R* and R spot = R^^/Rr/RZ. 
(Note, that for the isotropically emitting star, the stellar radius is 
related to the observed radius as R t = Roc yl — R B /R*, and is 
smaller than that.) 
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Figure 15. Constraints on the inclination and magnetic inclination of the 
system. The lines i = 30° and i = 85° are the lower and upper limits on 
the inclination of the system (M02). A small, p = 0, antipodal spot is not 
visible in the region cos(i + #) > k (below the corresponding line). A finite 
p shifts this line further to the left by angle p. The primary spot is eclipsed 
sometimes when cos(i + 6) < — n (to the right from the corresponding 
line), and, therefore, is forbidden. Upper hatched area gives the allowed 
regions of parameters i and 6, where the antipodal spot is blocked by the 
disc of inner radius R ln = 40 km, and the primary spot is visible all the 
time. Contours of constant variability amplitude A = 0.045 are shown 
for a black body spot (dashed curve, computed using Eq. |S) and for the 
spot radiation pattern corresponding to the 5th scattering order in a slab of 
Thomson optical depth rx = 1.5 (dotted curve, numerical simulations for 
a pulsar of frequency 435 Hz). Neutron star of mass M = 1.4M0 and 
radius i?* = 3i? s is assumed. 



8.3.3 Constraints from the absence of the secondary maximum 

The pulse profile of XTE J1751-305 is almost sinusoidal. This 
means that we see the emission from only one polar cap, as an 
antipodal spot would have created a secondary maximum (or a 
plateau) in the profile (e.g. B02, IViironen & Poutanenll2004 . At 
the same time, the primary spot should never be eclipsed by the 
star. The second shock can be obscured by the star or by the op- 
tically thick accretion disc. Let us consider first the possibility of 
obscuration by the star. 

Obscuration by the star. The condition of obscuration of the 
secondary spot (class I in classification of B02) is cos a s < 0. Us- 
ing analytical formula for light bending ^2} [ip — ip B is now the an- 
gle between direction to the observer and normal to the secondary 
spot, cos ip s = — (cos i cos 9 + sin i sin 9 cos <f>), and cj> is the pul- 
sar phase], we get cos(i + 0) > K, where k = R s / (R, — Rs). For 
a more compact star, _R* < 2R S (or about 1.8i? s if exact bending 
formula is used), such a region does not exist at all. 

In Flg.[B] we show an allowed region at the i — 9 plane where 
obscuration is possible. We restrict the inclination of the system to 
lie in the interval 30° < i < 85°, as argued by M02. The hatched 
lower region corresponds to the allowed space of parameters i — 9. 
One sees that the allowed region is rather small. It becomes even 
smaller, if we account for the finite spot size. The obscuration con- 



dition is then cos(i + 9 + p) > k, and the boundary of the region 
shifts further to the left by angle p (i.e. by ~ 15°), further dimin- 
ishing the allowed region. This region is smaller for more compact 
stars. Thus, we find that obscuration of the second spot by the star is 
highly improbable, unless the inclination is smaller than 30° and/or 
the stellar radius is larger than 3R S . In that case, the inner disc ra- 
dius has to be at least ~ 200 km, not to overproduce reflection. 

The primary spot is eclipsed sometimes in the upper right cor- 
ner of the diagram, cos(z + 9) < — n (i.e. i + 9 > 180° — arccos k 
- class III in B02 classification), and, therefore, is forbidden. Again 
for a spot of finite size, the limiting line shifts to the left by p. 

Obscuration by the disc. Let us consider now obscuration of 
the secondary spot by the accretion disc. First, we need to compute 
radii i?disc at which photon orbits (reaching the observer) cross the 
disc plane. One can easily show that in flat space 

R disc > R, Vl + cos 2 9 1 cos 2 i. (7) 

In Schwarzschild metrics, photon orbits cross the disc at somewhat 
larger radius due to the light bending. We computed numerically 
photon trajectories for every pulsar phase and found the minimum 
disc crossing radius. Now we can determine the region in the i — 6 
plane where this radius is larger than the disc inner radius which, 
as we argued above, is Ri n ~ 40 km (see Section l8.3.H . We show 
this region in Fig, ll5l (upper hatched region, upper right boundary 
is determined by the condition of the visibility of the primary spot). 
We see that second spot is not visible only for inclinations above 
~ 67° (for i?* = 3Rs). Larger R[ n will further reduce the allowed 
region. 

If we account for the finite size of the spot, the allowed region 
shifts to the left by the angular size of the spot. However, the lower 
limit on the inclinations is not affected much, because now it is 
a rather flat function of 9. We conclude that the second spot can 
be blocked by the accretion disc. We obtain constraints i > 67° 
(i > 70°) and 9 < 30° (9 < 35°) for R t = 3R S (2.5R S ) for the 
spot size of p ~ 15°. 

8.3.4 Constraints from the variability amplitude 

Important constraints on the physical parameters of the system are 
given by the variability amplitude. The observed peak-to-peak am- 
plitude A w 0.045 gives the minimum intrinsic variability. Let us 
first assume that all the emission comes the spot. For a slowly ro- 
tating star, the variability amplitude (for a black body spot) can be 
estimated from (B02; PG03): 

A = U/Q, (8) 

where U = (1 — R B /R r ) sini sin#. The dashed curve in Fig. 1151 
shows the relation between i and 9 satisfying the observed A. 

In reality, the star rotates rapidly and the radiation pattern from 
of the Comptonized emission is far from the black body. Rapid 
rotation causes a slight shift of the maximum emission from the 
phase where projected area has the maximum towards the phase 
where Doppler factor has a maximum (i.e. quarter of the period 
earlier, see PG03). Additionally, light travel time delays skew the 
profile. 

Change in the emission pattern causes even larger changes to 
the oscillation profile. The angular dependence of the radiation es- 
caping from a slab of Thomson optical depth tt = 1.5 (as sug- 
gested by spectral fitting, see Section 14.31 and Table |4j is shown 
in Fig. ll6l (see also Viiro nen & Poutanenll2004t) . It depends on the 
scattering order n which is in turn related to the photon energy as 
E/3kT sccd w (1 + 4kT c /m c c 2 ) n . Thus we expect 3-7 scatterings 
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Figure 16. Intensity (normalized to unity at fi = 1) of the radiation escap- 
ing from a slab of Thomson optical depth t"t = 1.5 for different scatter- 
ing orders n as a function of fi (here fi is the cosine of the zenith angle). 
The seed photons with intensity I(fJ.) = 1 are injected from the bottom 
of the slab. C alculation are done in Thomson scatt ering approximation (see 
ISunvaev & TitarchulJl 985 ; Viironen & Poutanen 20Q4|)- 



to contribute most of the flux in the PCA range. As an example we 
compute the light curves expec ted for n = 5 using fully rela tivistic 
method described in PG03 and Viiro nen & Poutanenl 120041) . Con- 
tour plots at the i — 8 plane for A — 0.045 are shown in Fig. 1151 
by dotted curves. We see that the Comptonized emission produces 
smaller variability compared to the black body, and for the same i, 
larger 6 is needed to get the same amplitude. 

If we take the upper hatched region as the allowed region 
of parameters, we immediately get that the magnetic pole is mis- 
aligned from the rotational pole by a very small angle, 9 — 3 — 4° . 
This angle becomes larger if it is only a fraction of the emission 
that comes from the spot. 

8.4 Distance 

We used the distance of 8.5 kpc throughout this paper, assum- 
ing that the source is close to the Galactic centre. M02 placed 
XTE J1751-305 at a distance of at least 7 kpc, using indirect ar- 
guments. There are, however, clues that this distance might be ac- 
tually much smaller. 

Firstly, the source is very bright, I/ pca k = 0.131/Edd, when 
placed at 8.5 kpc. It is about an order of magnitude brighter than 
SAX J1808.4— 3658, while all other properties are very similar. At 
higher accretion rates, typically more that a few per cent of the Ed- 
ding ton limit, atoll sources switch to a soft (banana) spectral state 
(e.g. lHasinger & van der Klis 1989). This is not the case here. But 
the transition luminosity can be quite high, in particular during the 
onset of the outburst, where hysteresis effect has been observed and 
the source can remain in th e hard state despite high accretion rate 
iMaccarone & CoppilEo Q3). For example, Aql X-l and 4U 1705- 
44 have been obse rved in the hard islan d state at lu minositi es ~ 2 
per cent of L E dd iBarret & Olivel2002tlD'one & Gierlihskill2003h . 



Therefore, having XTE J175 1-305 in the hard state at 13 per cent 
of the Eddington luminosity cannot be ruled out. 

ICumming. Zweibel & Bildstenl 1200 ll) argued that the crucial 
difference between the millisecond pulsars and non-pulsating atolls 
was in the average accretion rate (related to the size of the disc, 
which is smaller in the millisecond pulsars). In atolls, the accretion 
rate is high enough to bury the magnetic field under the surface. 
The diffusion time-scale of rebuilding the field is ~1000 years, 
much larger than a typical period between the outbursts. In mil- 
lisecond pulsars the accretion ra te, m, is never high enough to bury 
the field. ICumming et alj 1200 ll) estimated that magnetic screening 
is ineffective for m < 0.02mEdd- Then, 0.13 Lgdd would have 
been too much for the magnetic field to survive in XTE J1751- 
305 during the outburst. It would need D < 3 kpc for m to be 
low enough for the field not to be buried. This issue would require 
further study to resolve. 

If XTE J175 1-305 were located at D = 3 kpc indeed, it would 
decrease the blackbody apparent areas quoted in Table|4]by factor 
(8.5/3) 2 « 8 and disc inner radius by factor 8.5/3 « 3. The 
model DBPS (and DBTH), where the seed photons were from the 
blackbody component, would give the blackbody area of ~ 100 km, 
and could not be rejected on the spot size basis. The inner spot area 
in DBPF model would be ~3 km 2 and its linear size derived in 
Section l5.3.2l would be less by factor ~3. 



9 SUMMARY 

We have analysed the 2002 outburst of the second accretion- 
powered millisecond pulsar XTE J1751-305. The broad-band 0.7- 
200 keV spectrum of XTE J175 1-305 obtained by XMM-Newton 
and RXTE is complex. The best-fitting models consist of two soft 
components and thermal Comptonization. Most likely the cooler 
soft component comes from the accretion disc, the hotter one from 
the surface of the neutron star, while Comptonization takes place in 
a shocked region in the accretion column. We estimate the electron 
temperature of the plasma 30-40 keV and Thomson optical depth 
t ~ 1.5 (for a slab geometry). In our best spectral models, the tem- 
perature of the seed soft photons for Comptonization is higher than 
the temperature of the visible soft blackbody component which is 
consistent with a picture where photons from a heated spot under 
the shock are scattered away by the shocked plasma. In this model, 
the 'observed' (at the infinity) area of the blackbody emission is 
about 100 km 2 corresponding to a spot of radius ~ 5-6 km. The 
area corresponding to a shock was estimated to be 20 km 2 , which 
translates to a radius of 2.5-4.5 km depending on the assumed com- 
pactness of the neutron star and the viewing angle. 

The analysis of the RXTE spectra obtained during the outburst 
shows that temperature of the soft blackbody photons as well as 
the temperature of the seed photons decreased monotonically. The 
spectral slope of the hard Comptonized component did not vary 
much. We argued that this favours a constant geometry (e.g. slab) 
emission region and two-phase models, where energy is dissipated 
in the hot phase and the seed photons for Comptonization are the 
result of reprocessing of the hard X-rays. 

The pulse profile cannot be fitted by a simple sinusoid, a sec- 
ond harmonic is required. The mean peak-to-peak amplitude of the 
first and second harmonic is about 4.5 and 0.15 per cent, respec- 
tively. There is a clear energy dependence of the profile, the ampli- 
tude gradually decreases and the peak at higher energies appears at 
earlier phases (soft lags). The time lags reach ~ 100 /is and seem 
to saturate at about 10 keV. This behaviour is almost identical to 
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SAX J1808.4-3658 where lags are factor of two larger. The ob- 
served energy dependence of the pulse profiles can be explained by 
a model where blackbody and Comptonization components vary 
sinusoidally with a small phase shift. Such a behaviour can result 
from the Doppler boosting and different angular distribution of the 
emission from these components (PG03). 

The non-detection of the reflected component in the time- 
averaged spectra and the absence of the emission from the second 
antipodal spot (we see almost sinusoidal variations) put constraints 
on the geometry of the system. We argued that the inner radius of 
the optically thick accreting disc is about 40 km. In that case, the 
secondary can be blocked by the accretion disc if the inclination of 
the system is larger than ~ 70° . Blockage of the secondary by the 
star itself is highly unlikely for the neutron star radii i?* < 3i? s . 
The observed variability amplitude constrain the magnetic pole to 
lie within 3^4-° of the rotational pole, if most of the observed X-ray 
emission comes from a hotspot and a shock. 
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